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Abstract 

We study the structure of shock-free solutions of the compressible 
Euler equations with large data. We describe conditions under which 
the Rarefactive/Compressive character of solutions changes, and con- 
ditions under which the vacuum is formed asymptotically. We present 
several new examples of shock-free solutions, which demonstrate a large 
variety of behaviors. 
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1 Introduction 

We consider hyperbolic systems of conservation laws in one space dimension, 

ut + f{u)., = ^, ueW", (1.1) 

It is well known that, due to the absence of dissipative effects, classical 
(C^) solutions cannot be sustained, and generically, gradients blow up in 
finite time. This is a physical effect which is manifested by the development 
of shock waves, at which the conserved variable u becomes discontinuous. 
Once a shock wave forms, one must study weak solutions, and the analysis 
becomes much more difficult. 
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Breakdown of classical solutions for scalar equations is classical, go- 
ing back to Bethe and Hopf [12j, and was resolved for 2x2 systems by 
Lax [14^ I15j . These results state that in the presence of genuine nonlin- 
earity, nontrivial small data leads to gradient blowup in finite time. For 
larger systems, similar results are available, again provided the initial data 
is small [HI [21]. For 2x2 systems, a pair of Riccati-type equations of 
the form w = w"^ which blow up in finite time can be derived. In [10], John 
derives an analogous system of equations for the gradient variables, also 
with quadratic inhomogeneous part. These quadratic terms represent inter- 
actions between different nonlinear fields, including self-interaction terms. 
For small data, after an initial period of nonlinear wave interaction, the so- 
lution is essentially decoupled into n waves, each propagating with its own 
wavespeed. Each of these waves can be approximately treated as scalar, 
and so breaks down in finite time. When the data is large, some results are 
available for particular systems [H [3l [5] , but in general the breakdown of 
solutions with large data remains an open problem. For results in higher 
dimensional settings, see [21 1 [25l [26] . 

In this paper, we study the 3x3 system of Euler equations of gas dy- 
namics, which has one linearly degenerate field. The equations, representing 
conservation of mass, momentum and energy, respectively, are 

Pt + {pu)x' = 0, 
{pu)t + ipu^ +p)x> =0, (1.2) 
(ipti^ + pe)t + pu^ + up)x> = 0. 

We use a Lagrangian frame, co-moving with the fluid, given by x = J pdx'. 
The equations become 

n-Ux = 0, 

Ut+Px = 0, (1-3) 
(iu^ + e)t + {up)x = 0, 

and these are equivalent to ()1.2p [6,!3l]. Here t = 1/p is the specific volume, 
p is pressure, u is fluid velocity, and e is the specific internal energy. By the 
Second Law of Thermodynamics, for classical solutions, the third (energy) 
equation can be replaced by the entropy equation, 

St = 0, (1.4) 

see [B]. The system is closed by specifying a constitutive law; for convenience, 
we consider a polytropic ideal 7-law gas. For classical solutions, in regions of 
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constant entropy, the first two equations of (jl.3p close, forming the p-system 
of isentropic gas dynamics |27] . 

In the first author independently derived a set of Riccati type equa- 
tions for the gradients of sound waves, used these to give a consistent defi- 
nition of the rarefactive and compressive character {R/C character) of the 
nonlinear sound waves, and gave conditions which guarantee shock forma- 
tion, for data of arbitrary size. By ()1.4p . the entropy is a linearly degenerate 
contact field, stationary for solutions. These equations are analogous to 
those of [ini [20] , and can be generalized to other physical systems |5] . 

The presence of a stationary entropy profile of moderate strength means 
that nonlinear interactions between fields occur on the same scale as self- 
interactions, and can lead to surprising behavior [28j. In particular, a sound 
wave can change its R/C character across a contact discontinuity (entropy 
jump). In [HIS], the entropy field is assumed to be continuous, and conditions 
which guarantee gradient blowup are given. In this paper, we study the 
growth of gradients of shock-free solutions with a varying entropy profile 
that can include jump discontinuities, and we consider data having large 
amplitudes. In particular, we will demonstrate the consistency of results 
of [1] for C'^ entropy profiles and [28] for contact discontinuities. 

For our calculations, it is convenient to consider shock-free solutions, 
in which the velocity u and pressure p are C^, while the entropy S{x) is 

except at finitely many points, and one-sided limits of S and Sx exist 
everywhere. If the initial data are shock-free, the solution remains shock-free 
until the first derivatives of u and p blow up [231 ^\ ■ The following theorem 
generalizes results of [28l S] . 

Theorem 1. If the entropy is non- decreasing, a forward R (resp. C) can 
change its character only if it crosses backward C waves (resp. R waves); a 
backward C (resp. R) can change only if it crosses a forward C (resp. R). 
Symmetric results hold if m{x) is non-increasing, provided the character of 
the opposite simple wave is reversed. 

We say that a solution is eventually noninteracting if all wave inter- 
actions end in finite time: that is, for t large enough, the solution consists of 
three regions, defined by outgoing backward, stationary and forward waves 
respectively. The forward and backward waves are either rarefactions or a 
single shock; however the profiles of the rarefactions and contact field need 
not be explicitly known. A Riemann solution is the most obvious example 
of a solution that is eventually noninteracting, but we will present several 
other examples. 
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We define an asymptotic vacuum as a characteristic with a vertical 
asymptote in the (x, t)-plane. Since the characteristic is determined by the 
equation ^ = ±c, we necessarily have p — )• along this characteristic. 
Hence, along this characteristic, the p, c and p all vanish as t — )• oo. Note 
that the uniqueness theorem for ODEs implies that the vacuum is not taken 
on in finite time unless it is present in the initial data, see [33\ [T6] . 

According to ^8], it is possible to carefully choose an oscillating entropy 
profile which supports shock-free, space- and time-periodic solutions to the 
Euler equations. Thus we do not expect to be able to prove definitive gra- 
dient blowup results if the entropy is non-monotonic, so we largely restrict 
our attention to monotonic entropy profiles. 

Theorem 2. Assume that the (variation of the) data is compactly supported 
and that the entropy is monotone. Then a globally defined shock-free solution 
is either eventually noninteracting or contains an asymptotic vacuum. 

In other words, in a monotone entropy field, if a solution continues to 
interact for arbitrarily long times and does not contain an asymptotic vac- 
uum (so in particular (jl.Sp fails), then a shock necessarily forms in finite 
time. 

In some cases, we can predict the appearance of an asymptotic vacuum 
by a condition on the initial data: an asymptotic vacuum can occur in the 
solution only if the Vacuum Condition holds, namely 

uo{oo) — uo(— oo) > m{—oo) zo(— oo) + m(oo) zo{oo). (1-5) 

Here m and z are canonical thermodynamic variables which are nonlinear 
transformations of the entropy and density, respectively, given in (j2.2p . ()2.3p 
below. We show that if the vacuum condition (jl.Sp holds, then there are no 
eventually noninteracting shock- free solutions. Note that this vacuum condi- 
tion is identical to the condition that predicts the existence of the vacuum in 
the solution of the Riemann problem [27] , and is asymptotically equivalent 
to the condition for an embedded vacuum in general initial data |33j . 

We then provide several examples of shock-free solutions. In particular, 
if the entropy profile is piecewise constant and increasing and the data is rar- 
efactive, we show that the boundary of the interaction region is necessarily 
a forward characteristic. 

We also analyze a case in which the entropy is nonmonotonic. The 
entropy consists of two contacts of equal strength. By restricting data to 
be between the contacts, we show that shock-free solutions are everywhere 
rarefactive. Moreover, in this case we analyze the long-time behavior of the 
solution. 
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Theorem 3. There are three possible long-time behaviors: asymptotic vac- 
uum; infinitely reflected waves which converge to the vacuum state; and 
infinitely reflected waves with non-vanishing wavespeed and density. The 
long-time behavior is determined explicitly by a bifurcation parameter ( de- 
termined by the initial data. 

In particular, we obtain interacting solutions which asymptotically ap- 
proach vacuum at the specific rate 

s-r = 0(l)(l + t)-(T-i)/(T+i), or r = 0(l)(l+t)2/(i+T), 

where s and r are the Riemann invariants. As a consequence, we note 
that the presence of an asymptotic vacuum is a stronger condition than 
development of vacuum as t — >• oo, and in particular, p.Sp can hold even 
though no asymptotic vacuum is present. 

We contrast the cases of monotonic and nonmonotoic entropies: when 
the entropy is monotonic, shock-free solutions must be eventually noninter- 
acting or contain an asymptotic vacuum, while if the entropy is nonmono- 
tonic, waves can be reflected between the contacts infinitely often. 

Finally, we briefly discuss solutions containing a single shock. We treat 
this as a free boundary problem with specific conditions on either side of 
the shock, and discuss global solutions of this type. 

2 Equations and Wave Curves 

We restrict our attention to a polytropic ideal gas, with equation of state 

e = CyT = — — and pr = RT, 
7-1 

so that 

p = ife^/^"T-^. (2.1) 

Here S is the entropy, T is the temperature, R, K, Cy are positive constants, 
and 7 > 1 is the adiabatic gas constant, c.f. The Lagrangian sound 
speed is given by 

We use the coordinates of [28j: that is, we define variables m and z so 
that 

c = mz'^, p = rnc dz and t = — — dz, (2-2) 
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for some d. In fact, it suffices to take d = and set 

m = CmC^^^^- and z = C^t'^^-^^I'^ 

with constants 

= (d - 1) 1^ and Crr, = ^/K^C-'^, 

and it is easy to cfieck tliat (12. 2p becomes 

. m'^z^+'^ , z^-'^ , , 

c = mz , p = — ; — - — , and r = -. (2.3) 

We sliall continue to refer to m as tlie entropy variable. 

In these coordinates, for solutions, p.3p are equivalent to 

c 

zt H = 0, 

m 

P 

ut + mczx + 2 — m^: = 0, (2-4) 
m 

nit = 0, 

the last equation being ()1.4p replacing the energy equation, valid for smooth 
solutions. 

2.1 Hugoniot curves 

We describe the shock waves (and contact discontinuties) by the Rankine- 
Hugoniot conditions, 

e[T] = -M, 

e M = [p], (2.5) 

e [iu2 + e] = [up], 

where is the shock speed, and the brackets denote the change across the 
discontinuity, as usual. Using the identity [ab] = a [b] + [a] b, where q = 
{qo + 5i)/2, the third equation of (|2.5p simplifies to 

e(N+p[T]) = 0. (2.6) 

Shock waves correspond to solutions of (j2.5p with ^ 7^ 0: to describe 
these fully, we solve ()2.6p and use (12. 5 p to determine ^ and [u]. Using 

pr , d+l 

e = and 7 = - — -, 

7 — 1 a — 1 
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equation ()2.6p becomes 



{d+l)p[T] + {d-l)[p]T = 0. 

We now use (|2.3p and solve to get 



mo / a — ^ Zo 



which we write as 



having defined 



^ = /(% (2.7) 

mo Zo 



fiZ)-=\—r-. for d—<Z<d^. 

We now use (|2.5|) to describe u and ^, namely 

[n] = ±7b]M and e = ±Vb]/M- 
Substituting and simplifying as above, we obtain 

ui — uq = ibmo zq ^(Z) and ^ = itmo Zq /i(.Z'), (2. 
where g and /i are respectively defined by 



■=\ldTTi i-z^-^ ■ 

Here we have labelled the states on opposite sides of the shock with the 
subscripts and 1, without explicitly referring to the left and right states, 
and these can be interchanged in this description. We choose the signs in 
(|2.8p and the value of Z by referring to Lax's entropy condition: that is, 
we require that the (absolute) wavespeed be larger behind the shock, which 
implies that Z > 1 if zi is behind the shock, and Z < 1 if zi is ahead of the 
shock; if zi is the right state, say, we require Z > 1 for a backward shock and 
Z < 1 for a forward shock. In this way we can describe both forward and 
backward shocks with similar equations, as in [34j. For a further analysis of 
the structure of the wave curves and the functions /, g and h, see [3]. 
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2.2 Stationary Solutions 

The contact discontinuities provide another class of weak solutions, being 
solutions of ()2.5p with wavespeed ^ = 0. It is clear that contacts should 
satisfy [u] = [p] = 0, while the entropy variable m can vary: in this case, by 
(j2.3p we have 

ml zl^"^ = m\ z\+'^, (2.9) 

which in turn determines the jump in z. It is well-known that the entropy 
is linearly degenerate and these waves are contacts which are stationary in 
a Lagrangian frame |27) . More generally, we can obtain stationary waves by 
allowing m (or z) to vary while fixing u and p as constants: this is easily seen 
by direct substitution into equations (|1.3p . That is, any time- independent 
states (z, ti, m) given in some region by 

m = m{x), u = U, m{xf z{xy+'^ = P, (2.10) 

with U and P constants, form a stationary wave solution to ()1.3p . 



2.3 Isentropic Flow 

A simpler 2x2 system, known as isentropic flow, is obtained when the 
entropy S (or m) is taken to be identically constant, and the third (energy) 
equation of (II. Sp is dropped, to give 

Tt-Ux = 0, 

ut+p^ = Q, (2.11) 

with p = p{t), also known as the p-system. 

The p-system is considerably simpler than the full Euler system because 
the equations weakly decouple. Indeed, for solutions, ()2.4p becomes 

c 

zt H Ux = 0, 

m 

ut + mcZx = Q, (2.12) 
so that the Riemann invariants, given by 

r = u — mz and s = u + mz, (2.13) 

respectively, satisfy 

= and s' = 0. (2.14) 

Here ^ and ' denote differentiation along backward and forward characteris- 
tics, respectively, 

' = dt- cdx and ' = dt + cd^. (2.15) 
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3 R/C Character of Solutions 



We briefly recall results from the authors' previous papers [281 H] describing 
the local rarefactive and compressive nature of solutions. 



In a constant entropy field For isentropic flow (j2.1ip . the Riemann 
invariants r and s are constant along characteristics, and simple waves are 
described by 

mr = mi, Ur - Ui = mi{Za - Zb), (3.1) 

where the subscripts /, r, a, b denote the states to the left, right, ahead of 
and behind the wave, respectively. Noninteracting simple waves are classi- 
fied as rarefactive or compressive according to whether the characteristics 
diverge or converge, respectively. For isentropic flow, this is determined by 
the profile of the Riemann invariants. Since in the full system, entropy is 
stationary before shock formation, the following results extend immediately 
to isentropic regions in full 3x3 flows: 

Definition 3.1. f2^] In a constant entropy field, the local R/C character of 
a solution is: 

Forward R iff st < 0, 
Forward C iff St > 0, 
Backward R iff rt > 0, 
Backward C iff rt < 0. 

Lemma 3.2. 1341 In a constant entropy field, if an interacting solution is 
C"^ , the R/C character of each wave is preserved along characteristics. 

Lemma 3.3. \lJi\ \35^ Assume the initial data ro(x) and sq{x) ofr and s are 
C^, and the initial entropy is constant. If —st or rt is negative somewhere 
in the initial data, then \ux\ and/or \px\ blow up infinite time. 

In [I4j , Lemma 13.31 relies on an a priori assumption that the solution 
stays away from vacuum; this assumption is removed in [35]. Thus gradi- 
ents will blow up (shocks will form) in finite time if and only if there are 
compressive waves in the data; see also [I6] . 



At a contact discontinuity Simple waves preserve their character in 
isentropic regions, as the (derivatives of) Riemann invariants propagate with 
the wave. However, waves may change type when crossing a contact discon- 
tinuity which separates different isentropic regions. According to ()2.9p . the 
change in variables is 

Ur = ui and Zr = nii zi Q, (3.2) 
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where we have set 

Q=f!!}L)iTT^ SO also ^=g^. (3.3) 

m zi 

with corresponding changes in the derivatives of Riemann invariants rt and 
St by (j2.13p . It follows that if forward and backward simple waves cross 
the jump simultaneously, then one of the waves could change character. 
Following [3] , we will call the contact discontinuity a 1-contact if the entropy 
decreases, mi > rur (so Q < 1), and a 3-contact if mi < m,. (Q > 1). 

Lemma 3.4. [23^ A nonlinear wave changes its R/C value at a contact 
discontinuity when and only when one of the following inequalities hold: 

Rin '^mit ^ff Qmk <ui <mizi, 

^in ^out iff i^l^l <ui < Qmizi, 

Rfn C'^t iff -Qmk <ui < -mizi, 

^tn ^tut iff -mik <ui < -Qmizi, 

where y ■= yt denotes the time derivative, the subscripts denote incoming 
and outgoing waves (or the side of the jump), and the superscripts indicate 
the direction of the wave: — is backward, + is forward. 

Note that the conditions of the lemma are mutually exclusive, so only 
one wave can change its character at any time. Moreover, for a fixed jump, 
a change in one wave is possible only if the opposite wave has the right 
character. 

Corollary 3.5. At a 3-contact (Q > 1), the backward wave can change from 
R to C ( resp. C to R) only if both the incoming and outgoing forward waves 
are R (resp. C); the forward wave can change from C to R (resp. R to C) 
only if both backward waves are R (resp. C ). Similar conclusions hold for a 
1-contact, but the character of the incoming opposite wave changes. 

For later use, we record the change of Riemann invariants across the 
jump: it follows easily from ()2.13p . ()3.2p that 

For non-isentropic smooth solutions In |4j, the first author provides 
an appropriate definition of the R/C character for the full (non-isentropic) 
Euler equations. Recalling (|2.15p . define the quantities 

a := —p' l(? = Ux + mzx + -^j^ mx z and 

■- -p' /c^ = Ux - mzx - -AifUxZ] (3.5) 
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these are multiples of derivatives of Riemann invariants. 



Definition 3.6. f^f^ The local R/C character in a solution is 

Forward 
Forward 



R 




a> 0, 


C 




a <0, 


R 




/3>0, 


C 




/3 < 0. 



For solutions, it is easy to show using ()2.3p . (j2.4p and (j2.13p . that 

st + ca = 0, rt-c/3 = 0, 

so this definition agrees with and extends Definition 13.11 

The following Riccati type equations describe the growth of gradients: 



Lemma 3.7. If the solution of /ll.3\) is C^, th 



en 



' ' ~ ki (A:2(3a + /3) + a/3 — a^) and 



a 



where 



Moreover, 



P' = ki{-k2{a + 3P) + a^-^^), (3.6) 



7 + 1^ 7-1 .OVA 

fci := — -ZT-i ana ko := —, zzm^. (3.7 

2(7-1) 7(7 + 1) ^ ' 



\a\ or |/3| oo iff \ux\ or oo. (3.^ 



We note that these equations are similar to those derived by F. John in 
[lOj . but were independently derived from a different point of view by the 
first author in Condition (j3.8p coincides exactly with formation of a 
shock wave. 



3.1 Global R/C Structure 

In a constant entropy field, the R/C character of waves is preserved, but in 
a varying entropy field, it may change. We describe conditions under which 
the R/C character of an interacting wave changes. Essentially, the only way 
a wave can change is if it is nonlinearly superposed with reflections from 
the interaction of opposite waves with the background entropy field; this is 
consistent with the changes across a contact described above. 
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Theorem 1. Suppose a solution satisfies the Shock-free Condition. If the 
entropy m{x) (i.e. S{x)) is non- decreasing, a forward R (resp. C) can change 
its character only if it crosses backward C waves (resp. R waves); a backward 
C (resp. R) can change only if it crosses a forward C (resp. R). Symmetric 
results hold if m{x) is non-increasing, provided the character of the opposite 
simple wave is reversed. 

It follows from the proof that this statement includes changes of type 
from zero strength waves to C or R. 

Proof. First suppose the entropy is C^, and consider a forward rarefaction. 
We consider the evolution of a along the forward characteristic, propagating 
through a field of non-decreasing entropy. Also, suppose that /? > along 
this characteristic, so that our forward wave crosses no backward compres- 
sions. Let r denote the forward characteristic, parameterized by to < t- 
We prove by contradiction that a{t) > on F. Suppose not, and let 
be the first time for which a{t) = along F. Since /3 > along F, ki and k2 
are non- negative, and a{t) > for to < ^ < by (|3.6p we have 

a > 3kik2a — kia^ for to < t < t*. (3-9) 

Denote 



where the integral is along F. Using the integrating factor e -^'o 3'^i'^2rft^ (|3.9 



~ — ft 3kik2dt 17 1 ff 3kik2dt 

a(t) := a(t)e •'*(> and k+ := kie-'^o , 

the ii 

yields 

Dividing by and integrating along F, we get 



q' > — /c+d^. 



< / k+dt — 



a{t) Jto a{to)' 
Since d(to) > and a{t^) = 0, we must have 

lim / k^dt = -\-oo, 

which contradicts the Shock-free condition. We conclude that a > on F, 
with the lower bound 

r* 3kik2dt s 



1 + a(to) /t* Kdt 
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Now consider a point at which the entropy is not C^, (actuahy a contact, 
since entropy is stationary). By the Shock- free condition, one-sided hmits 
of m and nix exist and u and p are C^, so our R/C variables a and /? and 
characteristic T are defined up to x = with well-defined one-sided limits. 
If m{x) is continuous at x^, then a = and /3 = ^ are also continuous, 
so the above argument yields a(t*+) > 0, and we continue the characteristic 
forward in time. If the entropy has a jump at x*, then Corollary 13.41 applies, 
and the conclusion of the theorem follows. 

The proofs of other cases are entirely similar, and omitted. □ 

Boundary R/C character structure Finally we consider the R/C struc- 
ture at the edge of the support of the entropy profile. The main case of 
interest is that of no incoming wave and an outgoing rarefaction, corre- 
sponding to initial waves being compactly supported and no shocks outside 
the support of the entropy, respectively. 

Lemma 3.8. Suppose that the solution satisfies the Shock-free condition. A 
wave emerging from a region of varying entropy keeps its R/C character as 
long as there are no incoming waves of the other family. If m{x) is non- 
decreasing, a forward R (resp. C ) emerging to the right reflects a backward C 
(resp. R) hack into the region of varying entropy. A backward R (resp. C) 
emerging to the left of the varying entropy reflects a forward R (resp. C) 
into the region of varying entropy. Similar results hold if the entropy m{x) 
is non-increasing, but the character of the reflected wave is reversed. 

Proof. We consider only the first case: a forward rarefaction emerging from 
the right edge xi of the varying entropy, which is non-decreasing. All other 
cases are similar. 

First suppose there is a jump at xi: by Def. 13.11 and (I2.13p . the states 
to the right of the jump satisfy 

rr = iir — mr Zr = and Sr = + rUr Zr < 0, 

so that Zr < 0, where y := yt denotes the time derivative. Applying (j3.2p 
and using Q > 1, we see that 

ri = ill — mi zi = nir i.^ (1 — l/Q) < and 
si = nir Zr {1 + 1/Q) < 0, 

so that the waves on the left are a forward rarefaction and backward com- 
pression. 
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Now suppose the entropy is smooth (C^) and increasing up to xi, but 
constant for x > Xi. Then also, for x > xi, we have 

a > and /3 = 0, 

since there are no incoming backward waves and the emerging forward wave 
is a rarefaction. By continuity, for x near xi and x < xi we have a > 0, 
/3 0, and, by (j3.7p . /c2 > 0. Thus the forward wave is a rarefaction for x 
near xi, and moreover, in this neighborhood, 

< 0, so also /? < 0, 

so the backward characteristic reflected back into the varying entropy region 
is compressive. □ 

4 Shock-free solutions with a single contact 

Our main results refer to interacting solutions which interact only for finite 
times or contain asymptotic vacuums. When there is a single entropy jump, 
say at X = 0, then these are the only possibilities for shock-free solutions. 
By first treating a single contact discontinuity, we avoid issues of multiple 
reflections of waves considered in later sections. 

Consider the interaction of smooth isentropic waves at a contact dis- 
continuity separating constant entropy values mi and m^. We describe the 
states on either side of the contact by 

zi{t) = z{0-,t), Zr{t) = z{0+,t) and ui{t) = Urit) = u{0,t). 

We find the R/C character of the incoming and outgoing waves by using 
(j2.13p to calculate the Riemann invariants, and differentiating these in time. 

Lemma 4.1. Suppose the interacting solution contains no shocks. The solu- 
tion is eventually noninteracting if and only if the Vacuum Condition { j.5j) 
is not satisfied. If the Vacuum Condition holds, the solution contains an 
asymptotic vacuum. 

Proof. Suppose that interaction ends in finite time T ^ 1, and use subscripts 
g and d to denote the constant states on the left and right of the contact 
after the interaction has completed, respectively, so that Ug = ui{T), etc. 
Since there are no shocks, the outgoing waves leaving the interaction region 
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must be rarefactions. Using (|3.ip . we relate these states to the extreme 
states (subscripted by ±00) as 

Ug li — 00 — f^ — OO (-2— 00 -^9)5 

Uoo-Ud = mooiZoo- Zd). (4.1) 

We now use (|3.2p to relate the states across the jump, 
Ud = Ug and Zd = nig Zg Q. 
Now since nig = m_oo and nid = m-oo, we get 

lioo 00 ni—QQ -2— 00 nioQ Zqq — (1 ni—QQ Zg <^ 0, 
so that ()1.5p fails. 

Now suppose that the vacuum condition holds. Then the interaction per- 
sists for all time, and because there are no shocks, the outgoing waves must 
both be rarefactions. For t > 0, trace forward and backward characteristics 
back from the contact at x = to the initial time t = 0, to define functions 
X-{t) < and X-^-{t) > 0, respectively. Thus the forward rarefaction start- 
ing at X-{t) meets the contact at t, etc., and since characteristics cannot 
intersect, we have x_ < and i;+ > 0. Now, since Riemann invariants are 
preserved on characteristics, we have 

ui{t) +niizi{t) = uo{x^{t)) + nil zo{x-{t)) and 

Ur{t) — nir Zr{t) = Uo{x+(t)) — nir ZQ{x+{t)), 

where {zq{x),uo{x)) is the initial data. Using (j3.2p . we conclude that 

no(x_(t)) - uoix+{t)) + mi zo{x^{t)) + nij. zo{x+{t)) 
= {1 + Q)niizi{t)>0. 

Now since (jl.Sp holds, at least one of x_(t) or x^{t) must converge to some 
finite as t — )• 00. It follows that the incoming characteristic beginning 
at X* (and also those starting further out) cannot meet the contact in finite 
time, so remains bounded for all time. Thus this characteristic has a vertical 
asymptote, and the solution contains an asymptotic vacuum. □ 

Corollary 4.2. If a shock-free solution is eventually noninteracting, then 
the vacuum condition M.5\) fails, whatever the entropy profile. 
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Proof. The proof proceeds as above and (j4.ip continues to hold, provided 
subscripts g and d refer to the states on either side of the varying entropy. 
For large times t > T, the solution restricted to the entropy profile is a 
stationary entropy wave, across which velocity u and pressure p are constant. 
Thus Ug = Ud and (j4.ip yields 

SO the vacuum condition fails. □ 

In isentropic flow, the asymptotic vacuum is produced only by the inter- 
action of two opposite rarefaction waves [191 [33]. When a contact is present, 
a single strong rarefaction can produce an asymptotic vacuum, since the 
leading edge of the rarefaction crossing the contact may produce a reflected 
rarefaction, and the interaction of the initial strong rarefaction with the 
reflected rarefaction can lead to the vacuum. According to Corollary 13.51 
one of the incoming waves must be a pure rarefaction, while the other can 
contain compressive regions which are changed by the interaction at the 
contact, or there could be no opposite incoming wave. 

To be specific, consider a backward rarefaction (initially compactly sup- 
ported in (0,oo)) interacting with a 3-contact (increasing entropy jump), 
rrir > mi, at x = 0, with no incoming forward wave. The profile of the 
outgoing waves is determined by the traces of the states on either side of 
the contact discontinuity. By ()3.2p and (13. ip . the initial data consists of 
constant states {zi,ui,mi) and {zm,ui,mr) satisfying 

rrir Zm = mi zi Q, 

together with a backward rarefaction given by 

^0(2;) - Ui = mr Zm - mr Zq{x), 

for a decreasing function zq{x) with zq{0) = Zm- If there is some such 
that 

tio(x*) - mr zo{x^) = ui + mi zi, 

that is, 

mrZo{x^) = — - — mizi, or zo{x^) = z^, 

then all backward characteristics beginning to the left of x^, cross the contact, 
while all those starting at x > x* asymptote, with corresponding states 
approaching vacuum. Note that the support of the vacuum in the limit 
t — 7- 00 is some interval of the form [0, x^]. 



16 



5 Shock-free solutions with Monotone Entropy 



We now prove our main theorem, which describes the structure of shock-free 
solutions with a monotone entropy profile. 

Theorem 2. Assume that the (variation of the) data is compactly supported 
and that the entropy is monotone. Then a globally defined shock-free solution 
is either eventually noninteracting or contains an asymptotic vacuum. 

Proof. Assume without loss of generality that the entropy is non-decreasing. 
We assume that the solution is shock-free, has no asymptotic vacuum, and 
is not eventually noninteracting, and derive a contradiction. Denote the 
interval on which the entropy varies by [xq, xi]. Since the data is compactly 
supported, there is some T 1 such that 

I3{xi+,t) = and a{xo—,t) = for every t>T. 

Next, since there are no shocks, by Lemmas 13.21 and 13. 3^ the outgoing waves 
must be rarefactions, so that 

a(xi+,t)>0 and (3{xo—,t)>0 for each t>0. 

Because there is no asymptotic vacuum, all forward and backward charac- 
teristics pass through the interval [xo,xi] of varying entropy in finite time. 
It follows from Lemma 13.81 that 

a{xi-,t)>0 and /3(xi-,t)<0 for t>T. 

For each t > T, denote the backward characteristic starting from the 
point {xi,t) by r_(t), and define 

= min{xo < x < xi\ /3{^+,t) < 0, > x, (^,r) S r_(t)}, 

so Xt:{t) is the first possible point on the characteristic that f3 becomes 
positive. The curve is continuous, since /3 is continuous away from 

contacts, and if /3 first changes sign at a contact, /3(xc+) < < /3(xc— ), 
then we have = Xc. In this case, continuity of follows since 

/3(xc— ) is a continuous function of a{xc+) and (3{xc+). 

By definition, we have /3 < on the right of the curve while also 

/3 = on the curve (except possibly at a contact). Away from a contact, 
we must have > 0, so that a < at by ()3.6p . At a contact. 

Corollary 13.51 also implies that the forward wave is C, so again a < 0. 
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Now, since the solution is not eventually noninteracting, there is some 
such that 

a(xi— ,t#)>0, and /3(xi — , < 0, 

by Lemma 13.81 We now trace the forward characteristic r+ back from 
(xi — until it first meets the curve x*(t) at the point say. We 

now have 

^(a^t'^t) — ^'^'^ a(xi — ,t#)>0, 

while also 

(3{x,t) < for < X < x^, {x,t) £ r+, 
which together contradict Theorem [TJ □ 

5.1 Piecewise Constant Entropy 

We can describe the structure of solutions in more detail if we make the 
further simplifying assumption that the entropy is piecewise constant and 
monotone non-decreasing. By this we mean that the entropy has finitely 
many jumps, while u and p remain C^. Our results apply directly to mono- 
tone non-increasing piecewise constant entropy with appropriate modifica- 
tion, and we expect that the extension to general shock-free solutions with 
varying monotone entropy holds with technical changes in the proofs. 

Lemma 5.1. Suppose the entropy is piecewise constant non-increasing. If 
the initial data are never forward compressive but somewhere backward com- 
pressive, then there are no shock-free solutions. 




Proof. If a backward C leaves the varying entropy field, then shocks form 
in finite time by Lemma 13.31 If not, because there are only finitely many 
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isentropic blocks, we can isolate the left-most isentropic block B containing 
some backward C. By Corollary I3.5| forward C can form only at a 3- 
contact when the crossing backward waves are C. Thus there is no forward 
C in block B or the blocks to the left of it. Hence the backward C in this 
block cannot be cancelled. Thus the backward C must be obstructed by a 
backward asymptotic vacuum in (isentropic) block B. But in this case, a 
singularity forms in finite time by Lemma 13.31 □ 

Lemma 5.2. Suppose the entropy is piecewise constant non-increasing. If 
the initial data are never compressive, then shock-free solutions are never 
compressive. For such solutions, we also have: 

(a) If the data is constant to the right of the entropy jumps, and if there is 
some initial forward R, then there are no shock-free solutions. 

(h) If the solution is shock-free and eventually noninteracting, then the upper 
boundary of the interaction region is a forward characteristic which does 
not end at an interior contact. 




X Xo Xl 



Figure 2: Shock- free eventually noninteracting solution with rarefaction data 

Proof. Recall from Corollary 13.51 that a forward wave can change from R to 
C at a 3-contact, only when the crossing backward waves are C. Thus, any 
backward C must appear earlier than forward C. Since there are no C in the 
initial data, by Lemma [5. H shock-free solutions are always non-compressive. 
This proves the first statement. 

We now note the following consequences of Corollary 13.51 

1. If there are no incoming backward waves while the incoming forward 
waves are rarefactive on some piece of the right-most 3-contact, then 
the reflected backward waves are compressive. If this happens, a shock 
necessarily forms in finite time, by above. 
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2. In a shock-free solution, a forward R cannot be cancelled since the 
backward waves are nowhere C. 



Statement (a) now follows from 1 and 2, because the initial forward R 
remains R until it meets the last contact, and so reflects a backward C. 

Finally assume the solution is shock-free and eventually noninteracting. 
Let T denote the maximum time at which waves cross the right-most entropy 
jump (at xi). We claim the forward characteristic traced back from (xi,T) 
is the upper bound of the interaction region. Since our shock-free solution is 
nowhere C, any backward R above this characteristic would reflect a forward 
R, which in turn reflects a backward C at xi. If this characteristic ended at 
an interior contact, a backward rarefaction would emerge, and later reflect 
another forward wave, a contradiction. □ 

Lemma 5.3. Suppose the entropy is piecewise constant non-increasing. If 
the initial data are nowhere compressive, then shock-free solutions contain 
an asymptotic vacuum if and only if il.5\) holds. 

Proof. By Corollary 14.21 and Theorem [2l we need only show that if there is 
an asymptotic vacuum, then (II. Sp holds. 

First, we claim that the right-most backward characteristic, denoted by 
X = ^{t) forms an asymptotic vacuum to the right of the entropy field, 
so that <I>(oo) > xi. If not, the interaction must end in finite time by 
Lemma 15.2( b). Next, for any time T, there exists a forward characteristic 
^ which crosses all 3-contacts and which meets xi at some > T, so that 



Moreover, by choosing T large enough, we can assume ^'(O) is to the left 
of the support of the initial data. We denote by ipQ, . . . the points at 
which ^ crosses the x-axis and contacts, respectively, as in Figure [3l Also 
set 00 = (^(0),0) and (t)i = 

Since there are no compressions, we have for any {x,t), 



and recall that n is C^. We write 

tto(oo) - Uo{-(X>) = u{(j)o) - u{(j)i) + u{(j)i) - u{lpn) + u{lpn) - u{lpo) 



^'(0) < xo and ^{t^) = xi, t^>T. 



2 



'X 



> 



> n(0o) - n(0i) + u{Tpn) - uiipo)- 



(5.1) 
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X Xq Xi 



Figure 3: Asymptotic vacuum 

Now we use (|3.ip to write 

m(0o) - ui4>i) = rrioo zo{oo) - ruoo z{(j)i), 
and, telescoping, we write 

u{lpn) - u{lpo) = u{'4)n) " H h u{lpi-) - u{lpo) 

= rUn z{^pn-l + ) - rUn z{lpn-) + • • • 

+ m-oo 2:(V'o) - in-oo z{ipi-) 
> - mnz{ipn-) + m-ooZo{-oo), 

where we have used 

ruk+i z{ipk+) = Qk ruk z{ipk-) > ruk z{ipk-), 

by ()3.2p . Here the are the intermediate entropy levels (with rriQ = m_oo), 
and Qk > 1 the corresponding entropy jumps. Equation ()5.ip thus becomes 

no(oo) - no(-oo) > moo zo{oo) + m_oo •zo(-oo) 

- moo z{(j)i) - rrin z{ipn-), 

and the last two terms vanish in the limit as — )• oo, yielding (jl.Sp . □ 

6 Examples of shock-free solutions 

When the entropy is monotonic, a nontrivial shock-free solution must con- 
tain an asymptotic vacuum or must be eventually non-interacting. Here we 
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describe a general method for constructing eventually noninteracting shock- 
free solutions. 

Begin by specifying a piecewise smooth (C^) compactly supported en- 
tropy profile. Fix some constant xq (say the location of an entropy jump), 
and specify data {z,u) or (r, s) on a compact i-interval at this xq. Treat 
this as Cauchy data and evolve it spatially in both forward and backward 
directions. The equations for spatial evolution form a 2 x 2 system with 
varying coefficients, while the entropy is smooth. At an entropy jump, the 
jump is resolved by the Hugoniot conditions (12. 9p . which for a 7- law gas is 
simply a 2 X 2 linear map (j3.2p . 

In order to obtain global existence, we require only an a priori esti- 
mate \n \ [23], That is, we obtain a shock- free solution as long as a and 
P remain bounded in the half-plane t > 0. Taking the trace of the solution 
on t = then yields non-trivial initial data which generates this nontrivial 
interacting solution. By finite speed of propagation, for any fixed x, the 
solution will be stationary for t large enough, so the solution is eventually 
noninteracting. 

Because of Lemma 15.2( b). all non-compressive shock- free solutions with 
piecewise constant monotone entropy profile can be generated in this way, 
and we expect that all such non-compressive shock-free solutions have this 
structure as long as the entropy profile is monotonic. 

Rarefactions with two monotonic contacts By way of example, we 
explicitly construct an eventually noninteracting solution consisting of rar- 
efactions with two increasing entropy jumps (3-contacts). For convenience 
we set m_oo = 1 and define the entropy profile by jumps Qq > 1 and Qi > 1 
at Xq and xi, respectively, see (|3.3p . 

We specify Cauchy data by r{t,xo—) and s{t,xo—) = 0, with 

r(t,Xo — ) := rt{t,XQ — ) > 0, supported on t £ [0,T]. 

This means that there is a simple backward rarefaction wave emerging from 
Xq—, and no incoming forward wave. Equivalently, by (j2.13p . we choose 

z{t,Xo-) = Z{t), u{t,xo-) = -Z{t) for t G [0,r], 

where Z{t) is a positive monotone decreasing function. 

By applying Corollary 13.51 at both jumps, it follows that all waves in 
the solution are rarefactions, so that s < and r > everywhere. In 
particular, the maximum and minimum values of z are taken on at (xq— , 0) 
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and (xi+,cxd), respectively. Since the solution is eventually noninteracting, 
by (j3.3p . it follows that we have the uniform global bounds 



< z{x,t) < Z*, (6.1) 

where the bounds are given by 

Z* > Z{0) and < Qf^ Z{T). (6.2) 
We recall Lax's estimate of gradient growth in the p-system [14| [T5] . 

Lemma 6.1. In a constant entropy field, the gradients of Riemann invari- 
ants of solutions satisfy 

where the integrations are along the forward and backward characteristics 
connecting the points A = {xA,tA) cLnd B = [xB,tB), respectively, and 
where 

For convenience we choose xi large enough (or T small enough) that the 
backward wave to the right of xq meets the entropy jump at xi in negative 
time t < 0, so that only the forward wave crosses the second jump in our 
region of interest. This is clearly possible because we have uniform bounds 
()6.ip for z, and hence for the wavespeed. For the same reason, the backward 
wave to the right of xi meets the t-axis in some bounded interval [xijX*], 
so that the support of the data for the corresponds initial value problem is 
[a;o,X*]. By the mean value theorem, integrating the characteristics from 
(xo+,r), we get 

xi — xq >ciT 

xi -xo = C2(r* -T), and (6.4) 

x^, - xi = C3 r^,, 

for some values c, where T^, is the time at which the forward characteristic 
from (xo+,r) meets xi. In fact, since the forward wave is simple, C2 is 
exactly c(z(T, xo+)). 
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It remains to show that the backward rarefactions focus in the half-plane 
t < 0. It is convenient to describe the Riemann invariants in terms of the 
data Z{t), as follows. First, note that 

r{xo-,t) = -s{xo-,t) = -2Z{t) > 0, 

so that, by ([331), 

r{xo+,t) = -{l + Qo) Z{t) and 
-sixo+,t) = -iQo-l)Z{t). (6.5) 

By construction, there are two cases to consider, namely, a backward 
characteristic from a;o+, and a forward characteristic from X()+ to xi— fol- 
lowed by a backward characteristic from xi+. For the first case, we apply 
([O]) with A = {xo+,t), and B = (x,0), with t < T and xq < x < xi. We 
require 

^ = ^i2(A,S)-|K(x,5)dx>0, 
which reduces to 

R{A,B) 
^ ' jK{x,B)dx' 

and which clearly follows if 

• 1 mini? 

Now consider a forward characteristic from ^ to C = (xi,ti), followed by a 
backward characteristic from C to D = (x,0). Since —s > 0, (|6.3p implies 
that 

and again applying (j3.4p . we get 

r(C+) = -s{C-) < -s{A+' ~ ^ 



< -zit) 



2 - ' '2i?(^,C7-) 
(Qi-l)(Qo-l) 



2i?(^,C-) 

As in the first case, the backward rarefaction focusses at t < provided 

RiC+,D) 



r(C+) < 



/ K{x,D)dx' 
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which certainly holds if 

Z(t) ^ (mini?)2 

^ (Qi - 1) (Qo - 1) (X.-xi) maxK' ^ 

It follows that if (j6.6p and (j6.7p hold, then our solution satisfies the 
required properties. We now further estimate (|6.6p and (|6.7p . We fix the 
bounds and Z* > .Z'(O) for the data Z{t), use (|6.4p to eliminate X* — xi, 
and simplify, to find constants Ki depending only on Z^, Z* and Qi, such 
that, if 

-Z(t)<min{-^, ^^|, (6.8) 



then ((62]) and ((62]) hold. 

Finally, we show consistency of the construction, as follows. First, fix 
Z* and Qi, and choose xq, xi and T such that 

ici — xo > Ci r, 

where Ci is an upper bound for the wave speed in the region (xo,xi). Now 
choose Z{t) such that (j6.8p holds and such that 

Z* > Z{0) > Z{T) > Q^ Q^^ Z,. 
It is evident that these conditions are consistent by further requiring, say. 



< -Z{t) < 



2 2 
7* /n d+l r\ d+1 7 



7 Non-monotonic Contact Discontinuities 

We have shown that if the entropy profile is monotonic, then shock- free 
solutions are eventually noninteracting or contain an asymptotic vacuum. 
Here, by example, we show that this need not be true when the entropy 
profile is non-monotonic. Our entropy profile is piecewise constant with two 
contacts. As in earlier sections, we characterize the jumps using (j3.3p . Thus, 
we place a 3-contact Qo > 1 at xq and a 1-contact Qi < 1 at xi; without loss 
of generality we assume Qq = Q = l/Qi- We assume also that there are no 
incoming waves on either side of the interaction region, so the interactions 
are confined to the strip xq < x < xi. 

We consider the initial value problem with data prescribed in the in- 
terval (xo, xi), and constant outside that interval. We note that if the data is 
anywhere compressive, shocks must form in finite time, as the compressions 
cannot be cancelled at the entropy jumps. 
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Lemma 7.1. Fix the entropy profile as described above, and assume the 
initial data is constant outside the interval [xq,xi]. Then the solution is 
globally shock-free if and only if the data is nowhere compressive. 

Proof. It suffices to show that rarefactive data produces global shock-free 
solutions. We ffist solve the initial boundary value problem in the region 

Cl = {{x,t) I Xo < X < xi, t > 0}, 

with boundary data prescribed by the requirement that there are no in- 
coming waves. We then resolve the states across the entropy jumps and 
propagate the solution outwards as simple (rarefaction) waves. 
We obtain the boundary conditions by setting 

s{xQ—,t) = s{xo—,0) and r(xi+, t) = r(xi+, 0), 



and solving the Hugoniot conditions (|3.4p . After simplification, we write the 
boundary conditions as 

s{xo+,t) = ]^ _^ g r{xo+,t) + y^Tq ■5(^0-, 0), 

According to Corollary 13.51 the waves that are reflected back into the 
domain are always rarefactive, so we obtain global bounds on the derivatives 
of Riemann invariants. 

We now apply the existence theorem of |16| . Chap. 6, Thm. 3.1, which 
states that there is a unique global solution in which is nowhere com- 
pressive, provided a priori bounds are satisfied. We obtain explicit time- 
dependent bounds below, which suffice for application of the theorem. 

Finally, having solved the IB VP inside the domain Q, we again apply 
p.4p to obtain the Riemann invariants outside the domain, which are rar- 
efactions by construction. These thus propagate for all times without form- 
ing shocks. □ 

It is clear that these solutions are not eventually noninteracting, as there 
is always reflected rarefaction. It follows that either the solution forms an 
asymptotic vacuum, or waves continue to reflect back and forth between xq 
and xi for all times. 

We analyze a general rarefactive solution, as follows. Starting from the 
corner (xi, 0), we trace the reflected characteristic through the solution. Let 
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Figure 4: Infinite reflection of characteristics 



tfc denote the time of the A:-th intersection of this characteristic with the 
boundary, as in Figure [H and use the subscript to denote the corresponding 
interior state, so Z2k = z{xi—,t2k) and 2:2^+1 = z{xQ+,t2k+i)- We use (j3.ip 
to describe the states at either end of the backward and forward character- 
istics, respectively, by 

U2k - U2k+l = m {Z2k - Z2k+l), 
U2k - U2k-1 = m {Z2k^l - Z2k)- 

Since the waves outside Q are simple, we obtain 

'"2fc+2 -U2k = -Q {z2k+2 - Z2k), 
771 

U2k+1 - U2k-l = -Q {Z2k-l — Z2k+l), 



where we have used (j3.2p to express the outside states in terms of the interior 
states. Eliminating we get the same equation for even and odd n, 

{l + Q)Zn+l = 2QZn + {l-Q)Zn-l, (7.2) 

a linear diff'erence equation for Zn- We obtain a linear equation because the 
simple wave description ()3.ip and the jump conditions ()3.2p are linear in u 
and for m constant. 
Setting z„ = A", we get 

A2 - (1 + ?7) A + r? = (A - 1) (A - r?) = 0, 

where we have set 

,-1^.(0,1). 

It follows that the general solution of (j7.2p is 

..=^ie^+."^. (7.3) 

1 — 77 \ — T] 
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Theorem 3. For entropy profile as given above, and rarefactive initial data 
prescribed on (xo,xi), there are three possible long-time behaviors: asymp- 
totic vacuum; infinitely reflected waves which converge to the vacuum state; 
and infinitely reflected waves with non-vanishing wavespeed and density. The 
long-time behavior is determined by the expression C, = zi — rj zq, with a bi- 
furcation at Q = 

Proof. Since the data is rarefactive, zi < zq, so the second term of (|7.3p 
is positive and decreasing with n. If C > 0, then z„ is defined for ah n 
and approaches (^/{l — ij) >Oasn— )-oo. This imphes that the waves are 
reflected inflnitely often with uniformly bounded wavespeed. 

If = 0, then Zn = rf' zq, which clearly converges to the vacuum as 
n — >• cxD. Moreover, since Zn is defined for all n, the waves interact with 
the entropy jumps infinitely often, and the n-th characteristic has nonzero 
speed, so the solution does not contain an asymptotic vacuum. 

Finally, if C < 0, then for some A^, ()7.3p yields zat < 0, which contradicts 
the physical requirement that z > 0. We conclude that zat cannot be defined, 
which means that the iV-th characteristic never meets the boundary, and 
we therefore have an asymptotic vacuum. □ 

Corollary 7.2. Solutions with nonmonotic entropy may converge to vacuum 
as t —)• oo even though they contain no asymptotic vacuum. In particular, the 
condition of an asymptotic vacuum is stronger than the vacuum condition 

Proof. We have seen that C = yields a vacuum in the limit as t — > oo for 
all xq < X < xi. We show that in case C = 0) P-Sp holds as an equality. By 
()3.2p . at xi and xq, respectively, we have 

m 

uo{oo) = u{xi-,0+), m(oo) zo(oo) = — zo 

and 

u{xo-,ti) = u{xo+,ti), mzi = Qm{-oo) z{xQ-,ti), 
where zq = z{xi — ,0+) and zi = z{xQ+,ti). Now, by (jS.ip we also have 

u{xo-,ti) - uo(-oo) = m(-oo) (2:o(-oo) - z(xo-, ti)), 
u{xi-,0-\-) - u{xo+, ti) = m{zQ- zi). 

Now, m(— oo) = m(oo), and we calculate 

no(oo) - uo(-oo) = m(oo) {zq{-oo) - ZQ{xQ-,ti)) + m (zq - zi) 

m 

= m(oo) zq{—(X)) — -— Zi — mzi + Q m(oo) 20(00), 
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which yields 

uo{oo) — uo{—oo) — m(— oo) zo(— oo) — m(oo) zo{oo) 

TYi in 
= ~g ^1 -mzi + {Q - I) — zq. 

This vanishes since zi = rj zq = zq . □ 

Corollary 7.3. At the bifurcation point C = 0, the solution asymptotically 
approaches vacuum at the rate 

zo (l + iri-' - 1)) '''' <z^<zJl + (1 - v')" 

so, in particular, z{t) ^ 0{1) {1 + t)~^/'^ , which implies 

T = 0{l){l + tf-^l'^, or p = 0{\){\+tYl'^-^. 

Proof. Since z is monotonic along characteristics, we estimate the interac- 
tion times tn by 

Xl-XQ ^ ^ _ ^ Xl-XQ 
C{Zn) c{Zn+l) 

SO that 

n—l n 

m ^ ^ m " 

k=o k=l 

If ( > 0, then the terms z^*^ converge to k = (1 — riY/C,"^ ^ 0, so the times 
do not converge and tn ^ Kn for n large. 
When C = 0) we get 

n—l n 

mzi mzi 
which simplifies to 

j^-nd — I ^ tnTn Zq ^ _^ r?^"*^ — 1 

77~"<^ — 1 ~ xi — Xo ~ ^ ri~'^ — 1 ' 

which in turn yields 

1 + (1 - v'^) < ri-^'^ < 1 + (r/-'^ - 1), 

Xi — Xq Xi — Xq 

where cq = c{zo) = tuzq. Since z„ = 77" zq, the corollary follows. □ 
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8 Solutions with One Shock 



We briefly analyze the evolution of a single shock, as follows. The shock 
curve and states on either side satisfy the Rankine-Hugoniot relations ()2.5p . 
We analyze the interacting shock by treating it as a free boundary problem 
and imposing the conditions that the flow on either side of the shock is either 
isentropic or stationary. This yields a nontrivial interacting shock which is 
however not difficult to analyze. 

We parameterize the shock curve, S = (x{a),t{a)), together with Cauchy 
data on either side of S, subject to the Rankine-Hugoniot conditions. Ac- 
cording to (|2.7p and (j2.8p . we have 

— = a. 



zo 
nil 

mo 
Pi 

Po 

Ui - Mo 

t{a) - to 



/(«), 

±mz^ g{a), 
±fri'zh{a), 

x(a) 



i.l) 



mz g{a) 



da, 



where 



m 



mo 



! + /(«) 



{d + l)pi 



mn 



l/(d+l} 



are the averages, and subscripts refer to opposite sides of the shock; here 
the parameter satisfies 1 < a < d^^^'^~^\ 

Once we specify consistent Cauchy data, we can extend the solution to 
a neighborhood around the shock by solving the Cauchy problem locally 
on either side of the shock. Because we implicitly assume the Lax entropy 
conditions, the Cauchy problem is non-characteristic as long as the shock 
has nonzero strength, and so general existence theorems apply [17\ \W[ [7]. 
Moreover, if we establish a priori estimates for gradients, global existence 
follows. 

We specify the Cauchy data by assuming that the flow on one side of 
the shock is stationary, while the flow on the other side is isentropic. We do 
this by specifying that 



ui{a) = Ui, Pi{a) = Pi, and mo{a) = Mq, 
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so that 1 refers to the stationary solution and to the isentropic solution. 
This leaves a single free function, namely x{a), all other quantities being 
determined by (fHT^ . ([821) • 



Lemma 8.1. There are globally defined interacting solutions containing one 
shock which separates the plane into two regions, such that the solution is 
isentropic in one region and stationary in the other. 

Proof. Local existence follows because the Cauchy problem is not charac- 
teristic. Assume the shock is backward, with stationary solution behind the 
shock, so ui = Ur is constant, and the sign choices in (18. ip are all — . It is 
clear that the stationary solution behind the shock is globally defined. To 
verify existence in the isentropic region, we must estimate the derivatives of 
the Riemann invariants. 

To calculate the Riemann invariants, we differentiate the invariant along 
S, to get 

dr dx dt dx ( c{a) 



da ^ da ~^ ^ da ^ da \ ^(a)/ ' ^ ^ 
where we have used rt — crx = and ^ = — and similarly 

ds dx dt dx ( c(a) \ , , , 

-r = ^^-r^^t-r = ^^-r U + tK - (8-4) 



da da da da \ ^{a 
We now calculate from (18.11) and (12.31) that 



)\//2 

mi(a) = Mo/(a), = /(a)^ (a zo(a))'+\ and 



TT ( ^ M '{d+l)PiY/^'^'^ l + /(a) 
Ui - uo[a) = -Mo 



h{a). 



We use these to solve for zq^o) and no(a), and plug in to (|8.3p and (|8.4p to 
get 

da da 
rx = -r'^W and Sx = —S[a), (8.5) 
dx dx 

for explicit functions TZ and S. 

We now restrict a to a compact subinterval of (1, d^^^'^^^^), which implies 
uniform bounds on all thermodynamic variables, and so also TZ{a) and 5(a), 
and choose x{a) so that r^, and Sx are small enough on S that (|6.3p yields 
finite bounds for t > 0. 

The other case, of stationary solutions before the shock, is similar. Pro- 
ceeding as above, we obtain (18. 5|) (with slightly different TZ and S), and we 
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now must further restrict our data to ensure that Sx > 0, so the outgoing 
forward wave is a rarefaction, and we must ensure that the backward wave 
does not focus in the halfplane t > 0. □ 

It is interesting to ask whether the shock persists. Disappearance of the 
backward shock occurs in the hmit o — )• 1"*", so we must demonstrate that 
this Hmit cannot occur for finite {x{a),t{a)). It is weh known that in the 
Umit of vanishing wave strength, we have 

l-£i4 = 0(a-l) as a^l+, 

see |27j . The functions TZ{a) and 5(a) remain bounded away from in the 
Umit a — >• 1, so we obtain 

_ 0(1) da 
^ a — 1 dx^ 

and in particular, if rx remains finite, by (j8.3p . there exist u > and 6 > 
such that 

dx 

(a - 1) — > for 1 < a < 1 + 5. 
da 

Now, for < e < (5, we have 



x{l + 5) - x{l + e) = [ 
Ji 



dx , 
— da 
da 



> / da 

Ji+e a - 1 

, S 

= V log > oo, 

e 



as e —7- 0, so X —7- — oo as a — t- l"*". This implies that the shock strength can 
vanish only at infinity. 
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